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Abstract. We analyze the Fubini-Furlan-Rosetti sum rule in the framework of covariant baryon chiral 
perturbation theory to leading one-loop accuracy and including next-to-leading order polynomial contri- 
butions. We discuss the relation between the subtraction constants in the invariant amplitudes and certain 
low-energy constants employed in earlier chiral perturbation theory studies of threshold neutral pion pho- 
toproduction off nucleons. In particular, we consider the corrections to the sum rule due to the finite pion 
mass and show that below the threshold they agree well with determinations based on &xed-t dispersion 
relations. We also discuss the energy dependence of the electric dipole amplitude Eq+. 
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1 Introduction 

The Fubini-Furlan-Rosetti (FFR) sum rule was derived in 
the sixties utilizing the soft-pion techniques of current al- 
gebra pQ. It relates the nucleon anomalous magnetic mo- 
ment to an integral over the invariant amplitude A\ of 
pion photoproduction 

K ,, s = ^_ /^i m 4+.<v )t = 0)) (1 ) 

if one utilizes the Goldberger-Treiman relation gA m N — 
FirgirN, with mjv the nucleon mass, <?a the axial- vector 
coupling constant, F v the weak pion decay constant and 
<7„.jv the strong pion-nucleon coupling constant. Further- 
more, k v — n p — K n and k s — k p + K n are the nucleon 
isovector and isoscalar anomalous magnetic moment, re- 
spectively. The FFR sum rule is exact in the chiral limit of 
QCD and thus all quantities appearing in Eq. (JTJ are to be 
understood in the limit of vanishing light quark masses, 
m q = 0. The FFR sum rule has recently been reexam- 
ined in Ref. |2J. In that paper 1 , pion mass corrections to 
the sum rule were considered (in terms of a discrepancy 
function A]Sf(v,t)) and numerical evaluations based on a) 
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1 See also that paper for references to earlier work on the 
FFR sum rule. 



dispersion relations and b) input from heavy baryon chi- 
ral perturbation theory (HBCHPT) were presented. It was 
pointed out that in the strict framework of HBCHPT the 
nucleon pole positions are slightly moved, which leads e.g. 
to an incorrect curvature of the discrepancy function for 
energies below the threshold. A similar behavior due to the 
shift of pole or cut positions in the 1/m.N expansion was 
observed already in the discussion of the Compton cusp 
at the opening of the pion threshold [3], the spectral func- 
tions of the nucleon isovector form factors 4 or the scalar 
form factor of the nucleon 0. Note, however, that the 
kincmatical factors leading to the corresponding poles or 
cuts need not be expanded in HBCHPT as it is discussed 
e.g. in Ref. [HJ. Clearly, in a manifestly Lorentz-invariant 
formulation of baryon CHPT such problems do not arise, 
see e.g. jSEEl- The purpose of this paper is two- fold: We 
analyze the FFR sum rule in the framework of infrared 
regularization (IR) of baryon CHPT |SJ and demonstrate 
that the energy dependence of the discrepancy function 
is correctly given. Second, we also take a closer look at 
the pion mass corrections to the sum rule, which can only 
be systematically calculated in chiral perturbation theory, 
and the related threshold multipoles in pion photoproduc- 
tion. Our calculation includes all terms at third order in 
the chiral expansion and in addition, also the fourth or- 
der polynomial terms in the photoproduction amplitudes. 
Our aim is to show that within IR baryon CHPT one can 
describe pion photoproduction above and below thresh- 
old and that the dispersive representation can indeed be 
used to pin down certain low-energy constants (LECs), as 
suggested in 0. It is well-known that the chiral expansion 
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converges best in the unphysical region where all momenta 
can be very small. Consequently, in such regions LECs 
can be determined to a good precision if a correspond- 
ing dispersive representation is available. A more refined 
treatment including all fourth order terms and fits to the 
existing low-energy data from MAMI will be relegated to 
a future publication. 

The manuscript is organized as follows: In Sec. [21 we 
briefly recall the formalism of pion photoproduction and 
collect some results for the pion mass corrections to the 
FFR sum rule derived in [5] . We also present an alternative 
way of looking at these. Sec. 01 contains the results on the 
FFR sum rule, the discrepancy function and the related 
electric dipole amplitude Eq + as well as the slopes of the 
P-wave multipoles at threshold. We demonstrate that one 
can indeed determine LECs from the amplitudes in the 
unphysical region and end with a brief outlook. 



2 Formalism 

Consider pion photoproduction off the nucleon by real 
photons with k 2 = 

1 (k)+N{p 1 )^Tr a (q)+N(p 2 ) , (2) 

where pi (p 2 ) is the four-momentum of the incoming (out- 
going) nucleon (N) and a an isospin index (a = +, 0, — ). 
For the discussion of the FFR sum rule, only the isospin 
and + channels are of relevance, i.e. the physical channels 
7P — * ir°p and — > n n. The corresponding S-matrix is 
given in terms of four invariant functions Ai (i = 1, . . . , 4) 
that depend on two kinematical variables. Throughout, 
we utilize the notation of our earlier work [Hj and refer 
to that reference for a detailed discussion of the pertinent 
formalism. These invariant amplitudes can be calculated 
in baryon chiral perturbation theory and have the generic 
form (we do not display isospin quantum numbers and 
kinematical arguments) 

Ai = Af om + A l ° op + Af , (3) 

where the Born terms subsume the coupling to the charge 
and the magnetic moment of the nucleon (note often the 
alternative notation of "pole terms" is used for these con- 
tributions - sometimes even calculated employing the pseu- 
doscalar pion- nucleon coupling). All further counter terms 
are collected in the polynomial terms Af 1 . The non-trivial 
loop contributions (after renormalization of the single nu- 
cleon properties) are collected in the A l °° v . In Ref. [S], 
the Ai were calculated to third order (leading loop order) 
in relativistic baryon CHPT. In that formulation, a vio- 
lation of the power counting through the nucleon mass 
term is manifest. However, from the integral represen- 
tations given in it is straightforward to isolate the 
so-called infrared singular part 5 that contains the chi- 
ral long-distance physics and leads to a one-to-one cor- 
respondence between the expansion in loops and small 
momenta/pion masses. This can e.g. be achieved by the 



prescription given in which we also will employ. Sym- 
bolically, it reads f Q dx — > (J^° — f^°)dx = I + R, with 
I and R the infrared singular (irregular) and the regu- 
lar part, respectively. For the study of the FFR sum rule 
and related aspects, we are interested in the amplitudes at 
small energies and momentum transfer and thus use the 
variables^ = (s— it)/4mjv and = — (s+it— 2m^-)/4mjv, 
which are odd and even under crossing s <-» u (for further 
notation, see [EJ)- From the crossing properties of the Ai 
and their low-energy properties as detailed in %, , one de- 
rives the following representation for the polynomial pieces 
(note in particular the low-energy theorem for A\ j2] that 
forbids a constant term) 

Af = a\v 2 + a\v B + . . . , 
Af = a° + . . . , 
Af = a\v + . . . , 

Af = a° + ... . (4) 

At third order in the chiral expansion, only the leading 
term of A± contributes jS], whereas all other terms writ- 
ten down in Eq. start at 0(q A ). The mapping between 
these subtraction constants and the low-energy constants 
(LECs) used in the heavy baryon calculations |10II11II12| is 
given in the appendix. We remind the reader that although 
the contribution of Af formally starts as M 2 at threshold, 
in the chiral limit the expansion coefficients are singular 
leading to the famous contribution to the low-energy the- 
orem for the threshold value of Eo + at next-to-leading 
order in the pion mass expansion [SJ. We note that Af 
only feeds into the P-wave slope P 2 and into the D-wave 
at threshold in such a way that it cancels in the FFR sum 
rule at finite pion mass. Therefore, the subtraction con- 
stant a 2 has to be determined completely independently 
of the FFR sum rule, say by fitting to the slope of P 2 at 
threshold. In the following, we use the third order IR rep- 
resentation of the pertinent one-loop graphs (which con- 
tains an infinite series of corrections in the heavy 
baryon framework) but also use the fourth-order subtrac- 
tion constants displayed in Eq. (0J. Based on that repre- 
sentation, we attempt a simultaneous description of the 
^-dependence of the FFR discrepancy function A p (v, tthr) 
(as defined below), the energy dependence of the electric 
dipole amplitude Eq + for neutral pion production off pro- 
tons and the P-wave threshold slopes as extracted in ^Hl ■ 
Similarly precise information is not available for the neu- 
tron, therefore in the following we mostly concentrate on 
the proton. 

In Ref. |2 , the FFR sum rule was considered for finite 
mass pions and the following representation in terms of a 
discrepancy function An was derived: 

K N T 3 + A N (v, tthr) 

_ Am 2 N f°° di/_ v'lm A[ N >* a \ V \t = t thr ) 
eirg^N J UthT v V 2 - v 2 

A N (v,t = i t hr) 

= ^Hk t = t thl ) + Af(u, t = to*)) . (5) 
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A few comments on this equation are in order. First, the 
left-hand-side of the FFR sum rule gives the anomalous 

magnetic moment in the chiral limit, k n ~k n +0(ml /2 ). 
Therefore, if one uses the physical value of kn as in Ref . , 
one must include the corresponding loop and counter term 
corrections into the discrepancy function. However, for 
studying the pion mass corrections to the FFR sum rule, it 
is more appropriate to work with the chiral limit values of 
k p and K n , as discussed below. Second, in the chiral limit, 
t — > 0. This value can, of course, not be achieved in the 
physical world. We follow [2] and present our results at the 
minimal (threshold) value of t, t t hr = — M%/(l+M w /m) = 
—0.016 GeV~ 2 . Third, the Goldberger-Treiman relation is 
no longer exact at finite pion mass, to the order we are 
working, it takes the form |13U14| 
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with dis a LEC. As long as one only works at the physical 
pion mass, this effect is taken care of by utilizing the phys- 
ical value for the pion-nucleon coupling and the nucleon 
mass. If one, however, also wants to study the pion mass 
dependence of A at, as will be done here, one explicitly has 
to include this pion mass dependence. However, this effect 
only shows up in the terms cubic in the pion mass, which 
will not be considered in detail here (for a more detailed 
discussion of this topic, see e.g. |T£|). Note this pion mass 
dependence can only be systematically calculated in chiral 
perturbation theory and, eventually, in lattice QCD. 

We have explicitly worked out the quark mass expan- 
sion of the discrepancy function at threshold. For the pro- 
ton, it takes the form (modulo chiral logs) 



A p (v = Ufbr,t = tth r ) = a p M„ + PpMl + 



(7) 
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where £4 = mjv C4 and C4 = 3.4 GeV -1 |16j. The LECs 
eio5 and eio6 from the dimension four chiral pion-nucleon 
Lagrangian contribute to the proton and neutron anoma- 
lous magnetic moment at next-to-leading loop order (for 
a detailed discussion, see UJ). Their values are discussed 
below. Further, we have set the scale of dimensional regu- 
larization equal to the nucleon mass, A — m/y. Of course, 



Ap vanishes in the chiral limit. Notice the absence of chi- 
ral logs in the terms linear in the pion mass. The repre- 
sentation of Ap given Eq. (JJJ is exact to fourth order in 
the chiral expansion since it can be reconstructed from 
the HBCHPT results obtained in [TUIII1III2] . In addition, 
to arrive at these results, we had to include the small 
contribution from the slope of the D-wave combination 
D = M 2 + - E 2+ - P2- - E 2 -, that is D = D/q 2 at 
threshold. It is obtained from the invariant functions by 



D(s) 
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16 



dx (x 4 — 2x 2 + 1) .7-4(5, x) 



(9) 



with ^(s, x) a combination of ^2,3,4, see This contri- 
bution has not been calculated before. Note further that 

__ _ o 

Eq. J2JI includes the terms that renormalize n p to its phys- 
ical value, Kp, since we identify the left-hand-side of the 
FFR with the anomalous magnetic moment in the chiral 
limit. The pion mass expansion of A n (v = vthr,t = tthr) 
looks very similar, we find 
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(10) 



where the subtraction constants refer to the neutron am- 
plitude 771 — > ir°n as denoted by the superscript n. 



3 Results 

First, we must fix our input parameters. We use = 
92.4 MeV, M„+ = 139.57 MeV, M^a = 134.97 MeV, m p = 
938.27 MeV, m n = 939.57 MeV, g« N = 13.4, k p = 1.793, 
K n = —1.913. To the order we are working, the anoma- 
lous magnetic moment of the proton and the neutron are 
given in terms of two LECs from the dimension two chiral 
Lagrangian commonly denoted cq and C7 and loop correc- 
tions that start with terms of order ^7j. From that 
one can read off the chiral limit values for the proton and 
the neutron anomalous magnetic moments, 



Kp — 2.37 , K. n — 2. 



(11) 



in nuclear magnetons. For comparison, at third order in 

HBCHPT, we have A p = 2.85, A„= -2.98 QI|. 

The overall best fit to simultaneously describe the pro- 
ton discrepancy function, the electric dipole amplitude in 
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the threshold region and the three P-wave slopes is ob- 
tained with the following polynomial contribution to the 
invariant functions 

Af(v, v B ) = (191.3i/ B + 22CV 2 ) GeV~ 2 , 

Af{v,v B ) = -54.9 GeV" 4 , 

Af{v,v B ) = -155i/ GeV" 3 , 

Af(u, v B ) = 181.5 GeV" 3 . (12) 

While these coefficients appear large at first glance, in 
the appendix we show that they match quite nicely the 
corresponding LECs determined in fits to neutral pion 
photoproduction differential cross sections and the pho- 
ton asymmetry. Since the LECs can be understood to a 
good precision in terms of resonance saturation (excita- 
tion of the Z\(1232) and of vector mesons, see e.g. (10)1. 
the numbers appearing in Eq. l(T2"|l are indeed of natural 
size. The resulting threshold values for and the Pi 
(i = 1, 2, 3) are 

E 0+ = -1.19 [-1.23 ± 0.08 ± 0.03] , 
Pi = 9.67 [9.46 ± 0.05 ± 0.28] , 
p 2 = _9.6 [-9.5 ± 0.09 ± 0.28] , 
P 3 = 11.45 [11.32 ±0.11 ±0.34] , (13) 

in the conventional units of 10~ 3 /M 7T + and 10 _3 /M 2 + , re- 
spectively. The experimental numbers in the square brack- 
ets are from U^J- We note that Pi is obtained by adjusting 
the subtraction constant a^. The corresponding D-wave 
slope is 

D = 0.66 • 10~ 3 /M 3 + , (14) 

to be compared with 0.96 (0.92) from the MAID03 anal- 
ysis (the dispersive analysis of Ref. [T5]1- 

Let us look at these results in more detail. For the pion 
mass correction to the proton FFR sum rule, we find 

k p +A p {v = u thT , t = < thr ) = 2.18 (15) 

which should be compared with the left-hand-side of the 
sum rule, i.e. the proton anomalous magnetic moment in 
the chiral limit, cf. Eq. i|ll|) . Thus, within this framework, 
the FFR sum rule is fulfilled within 8%, which is of the 
expected size since pion mass effects are proportional to 
the small parameter [i = M^/jun — 1/7. Note again that 
this way of looking at the FFR differs from what was done 
in where the left- hand-side of the FFR was identified 
with the physical value of the proton anomalous magnetic 
moment and the discrepancy function defined there thus 

o 

differs from ours by terms ~k p —k p . 

Next, we display the resulting discrepancy function for 
the proton as a function of v at fixed t — £ t hr in Fig. ^ (for 
better comparison with Ref. |5j , this discrepancy function 
is taken with respect to the physical value of the proton 
anomalous magnetic moment). At v = 0, the discrepancy 
functions has indeed an extremum (as pointed out in jzj) 
and it increases with increasing v. We also find the pro- 
nounced cusp effect at v = f t h r , as it is expected. In con- 
trast to the prediction based on the MAID model (dashed 
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Fig. 1. The discrepancy function A v for the proton. Solid 
line: IR baryon CHPT. Dot-dashed line: third order relativistic 
CHPT 0. Dashed line: Dispersive result based on the MAID 
model, from [2|. Upper panel: v in the range from to 155 MeV. 
lower panel: Threshold region (135 MeV< v <150 MeV). The 
data points are calculated with the S- and P-wave multipoles 
from Ref. |18| (based on the formalism developed in |5j). 



curved), we do not observe any zero crossing for the val- 
ues of the subtraction constants given in Eq. I|12f) . Natu- 
rally, within our approach we should have a band rather 
than a line but we relegate a detailed error analysis to a 
later work, when we will also fit to all data of neutral pion 
photoproduction in the threshold region. Interestingly, the 
third order relativistic result from |S] (dot-dashed line) is 
very close to the IR result including the fourth order poly- 
nomial pieces. We also note that the incorrect behavior at 
v ~ observed in the HBCHPT calculation is of course 
not present in our manifestly covariant calculation, as it 
was suggested already in [2]. Note furthermore that we do 
not show Ap{y, i t hr) for values of v > 155 MeV because 
at v ~ 170 MeV, a steep rise due to the Z\(1232) reso- 
nance sets in In the lower panel of Fig. we focus on 
the threshold (cusp) region. We see that the relativistic 
and IR predictions are in good agreement with the MAID 
model and also with the data reconstructed form the mul- 
tipoles of Ref. ^Hj (for details, see again 0). 

We briefly discuss the pion mass dependence of A p<n ; 
for this purpose, we switch back to the full fourth order 
representation obtained from earlier heavy baryon results. 
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Consider the proton. We find for the coefficients in Eq. J2J 
a p = 2.19 GcV" 1 , 

[3 p = (-5.0 - 17.3) GeV~ 2 = -22.3 GeV~ 2 , (16) 

where we have used the physical values for and F n 
and the second term in [3 p is the contribution from the 
dimension four operators with the LECs eios and eioe- 
These counterterms determine the slope of k v and k s in 
the soft pion limit. Their values have been determined by 
using the fourth order formula for the anomalous magnetic 
moment from jHj as chiral extrapolation functions for the 
lattice QCD data of [52]. We get e W 5 ~ 0.45 , e W6 ~ 1.4 
at A = tun, which are of natural size. Since these are 
only very rough fits to the trend of the lattice data at 
too high pion masses, we refrain from assigning a theo- 
retical uncertainties to these numbers. Note that in case 
of the coefficient f3 p , we have large cancellations between 
the loop and the counterterm contributions proportional 
to the LECS a, which are -21.2 GeV~ 2 and 16.2 GeV~ 2 , 
respectively, for A = to at. We also stress that the first 
two terms in the quark mass expansion of A p at threshold 
give —0.13, which is to be compared with the full calcula- 
tion that gives —0.19, cf. Eq. (|15fl (where the contribution 
from the term ~ 2eio5 + e i06 is —0.34). For the neutron, 
we have no determination of the photoproduction counter 
terms and thus we can only give a n — a p = 2.19 GcV -1 
and $°°p = -18.23 GeV~ 2 . This value is comparable to 
the one of the proton. The contribution from the opera- 
tors ~ eios, eio6 is somewhat smaller than for the proton 
since for the neutron they appear with a different relative 
sign. 

Also well described is the energy dependence of the 
electric dipole amplitude Eq + (E 1 ) shown in Fig. [21 (with 
Ery = v— (t— M 2 )/4totv). This description of the data is as 
good as the complete fourth-order HBCHPT calculation 
(see e.g. Fig. 3 in This is not surprising since the 

third-order IR calculation generates all fourth-order HB 
corrections with fixed coefficients (from the expansion of 
the Dirac propagator) and these are the only fourth order 
loop contributions since loops with one insertion propor- 
tional to the dimension two LECs c$ do not contribute. 
Related to this is the observation that most of the non- 
trivial energy-dependence is given by the generalized cusp 
function [Tofe] 

E 0+ (E y )=a + b^jl-(^j , (17) 

where the parameter b is proportional to the charge ex- 
change scattering length a(TT°p — -> ir + n) (see also the de- 
tailed discussion in |2"Ull2"T] ). 

As noted before, similarly precise information on the 
neutron is not available. The predictions for the thresh- 
old multipoles are based on resonance saturation, the only 
experimental information is from the SAL experiment on 
•yd — > TT°d |23] that is consistent with the CHPT predic- 
tion |£o+% r | > l^o+ P thrl E2- We can nt to the energy- 
dependence of A™ as predicted by the MAID model and 
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Fig. 2. The real part of the electric dipole amplitude Eo+ in 
the threshold region. The solid line is the theoretical prediction 
as discussed in the text and the data are from Ref. |18|. 

the threshold multipoles as predicted by CHPT, but in 
the absence of more information on the energy dependence 
of e.g. the neutron electric dipole amplitude and any ex- 
perimental verification of the predictions for the P-wave 
slopes based on resonance saturation, the resulting num- 
bers are highly model-dependent and we thus refrain from 
showing them here. 



4 Summary and outlook 

In this paper, we have studied the FFR sum rule in the 
framework of covariant baryon chiral perturbation the- 
ory, extending some of the work presented in [2] . We have 
worked out the loop corrections to third order in the chi- 
ral expansion, corresponding to the leading loop contribu- 
tion, supplemented by the polynomial pieces up-to-and- 
including fourth order. Since the (sub-)threshold energies 
considered here are small, such a procedure is justified 
and would only lead to a readjustment of the subtrac- 
tion constants defined in Eq. if one includes also the 
fourth order loop graphs. We have shown that within 
this framework one can achieve a good description of the 
energy dependence of the discrepancy function for the 
proton, defined in Eq. J^J, together with the energy de- 
pendence of the electric dipole amplitude in the thresh- 
old region, cf. Fig. [2 and the P-wave slopes at thresh- 
old, see Eq. fH^ . The corresponding subtraction constants 
collected in Eq. (|12fl can be matched to the low-energy 
constants determined previously in HBCHPT studies of 
threshold neutral pion photoproduction and their result- 
ing values are of natural size (as detailed in the appendix). 
We find that the finite pion mass corrections to the FFR 
sum rule for the proton are small, of the order of 8 % 
(cf. Eq. |QSJ). As shown in Fig. ^ the ^-dependence of 
the discrepancy function for the proton has the proper 
behavior and agrees with the result of the MAID model. 
The unphysical behavior observed at v ~ in the heavy 
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baryon scheme [2] is absent in a covariant formulation 
as presented here. It is also interesting to note that the 
third order relativistic calculation free of counter terms 
from [E] also gives a good description of A p {v, t t in )- These 
findings corroborate the conjecture made in Ref. [2, that 
one can use the dispersive representation of the invari- 
ant amplitudes in the unphysical region to pin down low- 
energy constants of chiral perturbation theory (see also 
the appendix). It will be interesting to perform a com- 
plete fourth-order calculation and fit the corresponding 
LECs to the existing unpolarized and polarized thresh- 
old data of the reaction jp — > ir°p. This should further 
sharpen the conclusions made here. Work along such lines 
is under way |25|. 
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A Subtraction and low-energy constants 

Here, we give the mapping between the commonly em- 
ployed counter terms of the heavy baryon approach and 
the subtraction constants denned in Eq.(TJJ. To third or- 
der, one has only one P-wave counter term (with the LEC 
bp) that feeds into the multipole P3 



(18) 



At fourth order, there are two S-wave counter terms (which 
in the threshold region essentially act as one constant) 
|10|. The corresponding LECs a\ and are given by the 
following combinations of subtraction constants: 



127T a\ = 



m N 



34 
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Note that in the sum a\ + 0,2, the contribution from A2 
cancels (as noted earlier) and that these relations are to 
be taken at A = tun- Similarly, at fourth order there are 
two independent counter terms (with the LECs £1 and £2) 
that modify the P-waves Pi and P2 ^2 
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(20) 



At first glance, one might conclude from Eas. ll9l2UI) that 
there is a mismatch in the number of subtraction constants 



and counter terms. Note, however, that various subtrac- 
tion constants feed into the S- and the P-waves so that 
finally only two independent structures remain for the S- 
wave and two for the P-waves. 

It is interesting to compare the numbers derived from 
Eq. I|12(l with the earlier determinations of these counter 
terms in the HBCHPT framework. Note, however, that 
we did not include all fourth order loop corrections here, 
so that the values of the subtraction constants effectively 
subsume some of these effects. This is not the case for the 
LEC bp since it already appears at third order. Our value 
for dp translates into bp — 14.4 GeV -3 . This compares 
well with the third order fits of Ref. [HI], b P = (15.8 ± 
0.2) GeV -3 , and of Ref. |TT|, b P = 13.0 GeV" 3 . Note that 
the corresponding value in |12j comes out smaller due to 
additional loop effects. Consider next the LECs contribut- 
ing to the electric dipole amplitude in the threshold region. 
We get 47rmAr(ai + a 2 ) = 56.2 GeV -3 from the constants 
in Eq. O compared to 31.8 GeV" 3 , 77.8 GeV -3 and 
71.2 GeV -3 from (TH], Q31 and |T2], respectively. Individ- 
ually, we have at = -3.67 GeV" 3 and a 2 = 8.43 GeV" 3 , 
which is different from but comparable in size to the free 
and resonance fits to the various sets of Mainz and Saska- 
toon data (compare e.g. table 1 in 52)- As it was al- 
ready stressed in these earlier papers, the LECs a\ and 
C12 can not be well determined individually from fits to 
the data in the threshold region. Next, consider the P- 
wave P\. The determination of £1 = 16.6 in ^5] trans- 
lates into a\/2 — mjvaj — a% = 175.7 GeV -3 , which is siz- 
ably larger than the value of 59.8 obtained from Eq. \Yty . 
This is expected since in ^2] it was shown that there are 
large cancellations between fourth-order loop and counter 
term contributions, which we represent by the polynomial 
term only. This discrepancy is even more pronounced for 
the combination of subtraction constants that can be ob- 
tained from p2- Utilizing £1 = —19.7 from (T3| and the 
values from Eq. Ijl2|l . we obtain a\j2 + mjvaj — mjva!] = 
-208.5 GeV" 3 and -3.3 GeV -3 , respectively. This shows 
that the cancellations between the fourth order loop and 
counter term contributions are even stronger in P2 than 
in P x . 
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